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Abstract 

The degree distributions of many real world networks follow power-laws whose exponents tend to fall 
between two and three. Within the framework of the Barabasi- Albert model (BA model), we explain this 
empirical observation by a simple fact. To that end we propose a modified BA model with one parameter 
that serves as a regulatory factor for the growth rate of added links in scale-free networks. The regulatory 
factor has something to do with the obvious fact that one link has two nodes. The modified model also 
allows to connect nodes by newly added links that do not necessarily emanate from new nodes. Another 
related model using the master equation is also given, from which the same power-law degree distribution 
can be derived. 



1. Introduction 

In recent years, it has been revealed that many real 
networks obey the power-law degree distribution: 

p(k) oc fc~ 7 . 

Here p(k) is the probability that a randomly chosen 
node (or vertex) has k links (or edges), and the expo- 
nent 7 is a constant, usually 

2 < 7 < 3. 

Such networks are called scale-free in [1]. Some ex- 
amples of scale-free networks and their exponents are 
exhibited in the following table ([3]). 



Network 


7 


WWW (indegree) 


2.1 


WWW (outdegree) 


2.72 


Internet (domain) 


2.1-2.2 


Internet (router) 


2.48 


Movie actors 


2.3 


Phone-call 


2.1 



There exist also a few networks whose degree dis- 
tributions follow power-laws with exponents 7 < 2 
or 7 > 3. Traditionally, networks were thought of 
as random graphs, which were studied by Erdos and 
Renyi in the 1960s. For those the degree distribution 
is a binomial or Poisson distribution, which is peaked 



about the average and decaying exponentially away 
from it. 

In order to explain the above empirical fact of real 
networks, Barabasi and Albert proposed a model (BA 
model) in [1]. Their model could deduce the degree 
distribution: p{k) oc fc~ 3 , i.e., 7 = 3. 

Now we are able to find many other network mod- 
els by which power-law degree distributions with a 
variety of exponents can be deduced (see [4, 5] or 
references therein). However, it seems that the BA 
model even now provides a constructive model of com- 
plex networks for interested researchers and students. 

The objective of this paper is to produce scale- 
free networks with exponents between two and three 
based on the BA model. For that purpose we propose 
a modified BA model having one parameter. The 
parameter regulates the growth rate of added links in 
scale-free networks and it has something to do with 
the obvious fact that one link has two nodes. Another 
related model using the master equation ([5]) is also 
given, from which the same power-law degree distri- 
bution can be derived. 



2. A modified BA model 

In this section, reviewing the BA model ([1]) in 
more detail, we will propose a modification of the 
model. The modified BA model is involved with a 
parameter and it coincides with the BA model in a 
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particular case. We derive the condition for the pa- 
rameter and describe how it affects our model. 

The BA model is constructed by two ingredients. 
First, a network grows by continuously admitting new 
nodes to it and new nodes connect to nodes already 
present in the network. Thus it expands by acquiring 
more nodes and links over time. Second, the man- 
ner in which a new node connects to existing nodes is 
not completely random. The mechanism of connection 
is based on so-called 'preferential attachment'. This 
means that the probability that a new node connects 
to an existing node is proportional to its degree. 

The proposed model is as follows. We start with 
a network having a small number of nodes, too, and 
some links. Notice that the value too and the num- 
ber of the original links are not important, since we 
eventually take the limit with respect to time. New 
nodes having some new links constantly enter the net- 
work. Whenever a new node enters the network, we 
advance the time-step by one, and m new links are 
added to the network in each time-step, including the 
links emanating from the new node. Here to is a pos- 
itive integer. 

In our model, two nodes which already exist in 
the network are allowed to connect each other. This 
possibility was not considered in the BA model ([1, 
2]). Connections to existing nodes are always exe- 
cuted by preferential attachment. Preferential attach- 
ment means that the probability, V, that any node 
connects to a node v is proportional to the connectiv- 
ity k of that node v, so that 



of (2), we have 



V(k) 



(i) 



After t time-steps, we obtain a network with t + rriQ 
nodes and mt links plus the original links. Because of 
the mechanism of preferential attachment, a node that 
happens to acquire more links than others will con- 
tinue to increase its connectivity; a 'rich-get-richer' 
phenomenon. 

Thus the rate at which a node Vi with ki links ac- 
quires more links is given by 



dkj ^ ^ , , 



(2) 



where g is a proportionality constant, which is to be 
determined later. The denominator in (1), which we 
denote by 

C{t) = V hi 
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is the total of all degrees in the network at time-step 
t. In view of (1), taking the sums over i in both sides 



dC(t) 
dt 



9- 



(3) 



Since to links are added to the network in each 
time-step and each link is counted twice in counting 
degrees, the total degree C(t) at time-step t satisfies 



C(t) w 2mt. 



(4) 



Therefore, in view of (3) and (4), the growth rate 
g of the total degree is not equal to to as in [1, 2], but 
approximately equal to 2to. The growth rate will be 
equal to to, when one completely neglects the degrees 
of new nodes. On the other hand, it is equal to 2m, 
when one takes a full account of them. So it appears 
that taking g = to underestimates g. 

Considering the contribution of new nodes to the 
growth rate g to some extent, we introduce a param- 
eter (3 into the BA model by g = (3m. It is reasonable 
to assume that (3 satisfies 1 < (3 < 2 from the above 
argument. The parameter serves as a regulatory fac- 
tor for the growing network. 

Thus equation (2), together with (4) and g = (3m, 
can be written as 



dki ki ki 

~dt ~ m 2mt ~ 2t' 



(5) 



For the value of (3 we need the corresponding initial 
condition imposed upon a new node Vi, which enters 
the network at time-step U. This condition is the only 
equation that prescribes the new node, while (3) and 
(5) are supposed to be equations for those nodes that 
are already part of the network. 

From (3) we see that g = (3m is the amount of 
degree's increase for the part of the existing network 
within one time-step. On the other hand, 2to is the 
total amount of degree's increase. Thus the initial 
condition is given by 

h(U) = (2-0)m, (6) 

for each new node0 

Next let us examine the condition for (3 in further 
detail from (6). It necessarily implies (3 < 2 and (3 > 1. 
The former is obvious. For the latter, suppose (3 <1. 
Then it follows from (6) that the number of new links 
emanating from a new node is greater than to, which 
is impossible. Therefore, the appropriate condition for 
the parameter (3 is 

1 </3 < 2. 

Moreover, in the case of 1 < j3 < 2, notice that condi- 
tion (6) permits two nodes already present in the net- 
work to connect each other by some new links. The 



"(6) is the major correction to the previous versions. 
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number of such links, among to newly added ones, is 
equal to 

to- (2-/3)m = (J3- l)m. 

When /? = 1 as in [1, 2], the possibility is forbid- 
den. That is, in this case all new links are connections 
between only a new node and nodes of the existing 
network. 

In Section 3, under the condition 1 < (3 < 2, we 
show that this modified BA model, (5) together with 
(6), generates a scale-free network with exponent 



7 : 



1 



2 



So we have 2 < 7 < 3 for our model. In Section 4, 
using a quite different approach as in [5], a similar 
power-law will be obtained by means of the master 
equation based on (5), too. 



3. Scale-free networks with exponents 
between 2 and 3 

Under the assumption 1 < (3 < 2, the solution of 
(5) with initial condition (6) is given by 



/,',-(/) = (2- ( j 



(3/2 



(t > U), 



where ti is the time when a node Vi enters the network. 
It is assumed that the time interval s = t i+1 — t j (for 
each i) is a constant, usually s = 1. 

First we verify condition (4) for the model. The 
total degree C(t) equals 



(2 - P)m 



(3/2 



t\ (i/2 

V 1 



(3/2 



where ti is the maximum time-step such that ti < t. 
The sum in the square brackets is approximately equal 
to 2(2 — ft)~ 1 t by the continuous approximation: 



dx. 



Hence we see that (4) is satisfied. 

In order to obtain the degree distribution for this 
model, we could adopt the method employed in [2]. 
However, we will use the following alternative ap- 
proach, because it is easier to follow and, moreover, 
it gives the same asymptotic behavior for the degree 
distribution as the method in [2] . 

For the probability p(k) that a randomly chosen 
node v has a connectivity k, it suffices to estimate the 



number of time-steps U's for nodes i^'s, which satisfy 
the inequalities 



k-l< h(t) = (2- (3)m ( - 



0/2 



< k, 



(7) 



or alternatively 

k- 1/2 < ki(t) < fc + 1/2. 
Let us take (7) for the estimation and rewrite it as 



t 



k 



(2 - f3)m 



-2/0 



<U<t 



k-l 



(2 - f3)m 



-21(3 



In the case of s = 1, the number of nodes satisfying 
(7) is equal to the difference of both sides of the above 
inequalities: 

t{(2 - f3)m} 2 ^ 3 {(k - I)- 2 '? - k-W}, 



or 



t{(2 - p)m} 2 ^{(l - l/k)- 2 '? - l}k 



-2/(3 



Making use of (1 - 1/fc)- 2 /' 3 w 1 + 2/3~ 1 /k for large 
fc, we see that the number of nodes satisfying (7) is 
approximated by 

2(3~ 1 t{(2~(3)m} 2 ^k- 1 - 2 ^ 3 . 

Since at time t there exist t+mo nodes, the probability 
p(k) is thus given by 



2/TH{(2- (3)m} 2 'P 



t + to 

In the limit t — > 00 we conclude 

p(k) ocfc-( 1+2 /«, 
i.e., 7 = 1 + 2//3, for all large k. 

4. Another related model 

In this section we assume 1 < (3 < 2. Related to 
the model discussed in Section 2, we propose another 
by means of the master equation as in [5]. Making 
use of (5), we show that it also leads us to the same 
power-law. 

Let p(k, ti, t) be the probability that a node Vi, en- 
tering the network at time-step ti, has a degree k at 
time t. The master equation for it is given by 

p(k,t i ,t+l)=P^^p(k-l,t i ,t) 

+ (l-p^)p(k,ti,t), (8) 



3 



because the degree of a node with k links increases by 
one at the rate f3k/2t from (5). 

We define the limiting probability 



p(k) = lim \ V p(k,ti,t) 



t^oo t 

for deriving the degree distribution. It follows from 
(8) that the recurrence relation for this p(k) is given 

by 

(t + l)p(k) = (3^tp{k - 1) + (l - tp(k), 



or 



p(k) 



k-l 



p(fe-l). 



fc + 2//? J 

Using well-known properties of the beta function and 
the gamma function, it can be written as 

_ r(fc)r(i + 
p{k) ~ T{k + \ + 2/p) p[l) 

= B(fc,l + 2//?)p(l) 

for all large k. The last asymptotic behavior follows 
from Stirling's formula. 
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